Superfluid density in gapless superconductor CeCoIn 5 
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Temperature dependence of the London penetration depth A measured in single crystals of 
CeCoIns is interpreted as caused by a strong pair-breaking scattering that makes the superconductiv- 
ity in this compound gapless. For a gapless d-wave superconductor, we derive A = Ao/y 1 — (T/T c ) 2 
caused by combined effect of magnetic and non-magnetic scattering in excellent agreement with the 
data in the full temperature range and with the gapless s-wave case of Abrikosov and Gor'kov. We 
also obtain the slope of the upper critical field at T c that compares well with the measured slope. 

PACS numbers: 74.20.-z, 74.20.Rp, 74.70.Tx 
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The heavy-fermion superconductor CeCoIn 5 is still un- 
der intensive scrutiny after its discovery in 2001 
This is clean (the mean-free path exceeds by much the 
coherence length £ a b w 80 A), nearly two-dimensional 
(with small 3D pockets) 0|, d-wave superconductor 
In the normal phase for T > T c — 2.3 K, this mate- 
rial is a strong paramagnet [f|. The main interest of 
the community has been focussed on low temperatures 
and high fields where the inhomogeneous Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) phase is suspected to exist. 
In this work we arc interested in zero-field superfluid den- 
sity and the c directed upper critical field H c2 near T c , the 
domain removed from complications of FFLO and para- 
magnetic constrains. Understanding the ground state 
properties is of utmost importance for definite predic- 
tions about existence of more complex phases such as 
FFLO. 

Single crystals of CeCoIns were grown from an In flux 
by combining stoichiometric amounts of Ce and Co with 
excess In in an alumina crucible and encapsulating the 
crucible under vacuum in a quartz tube. Details of the 
crystal growth are described elsewhere [H. The crystals 
used were 1 x 1 x 0.2 mm 3 and magnetic measurements 
showed practically no hysteresis. 

The magnetic penetration depth was measured with a 
tunnel-diode resonator sensitive to changes in A of about 
1 A. Details of the technique are described elsewhere 
In short, properly biased tunnel diode maintains 
sclf-rcsonating tank circuit on its resonant frequency 
uj ~ 14 MHz. A sample is inserted into the coil on a 
sapphire rod. The changes in the effective inductance 
cause a shift in to proportional to the real part of the dy- 
namic magnetic susceptibility. The system is calibrated 
by matching the T dependent skin depth just above T c . 
To probe X a b, a small ac magnetic field (~ 20mOe) is ap- 
plied along the c-axis, so that the screening currents are 
flowing in the ab plane 

The London penetration depth X a b(T) as a function of 
temperature shown in Fig.[T] With excellent accuracy, 
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FIG. 1: (Color online) The data on A(T) and the fit to Eq. © 
with A(0) as the single fitting parameter. Grey line: data 
from Ref. Inset: low temperature part of A(T) replotted 
vs. T 3 / 2 and T 2 for comparison. 



the data from 0.5 K all the way to 2.1 K are described by 
Xab = A(0)/ Vl - t 2 , t = T/T c . (1) 

In fact, our data are close to those reported in Ref. 
(the grey band in Fig.[T]); our interpretation, however, is 
quite different. 

We note that the superfluid density 1/A 2 deviates from 
the zero-T value approximately as T 2 at low T's, un- 
like both the s-wave exponentially weak behavior and the 
clean d-wave linear T dependence. On the other hand, 
we recall that Abrikosov and Gor'kov (AG) in the semi- 
nal paper Q on pair-breaking by magnetic impurities in 
isotropic s-wave materials found the dependence |T]) in 
the whole temperature range from to T c for a strong 
spin-flip scattering when T c is suppressed to nearly zero 
and the superconductivity is gapless. 

This suggests that a similar situation takes place in 
CeCoIns, although the reason for the pair-breaking may 
not be the spin-flip scattering on independent magnetic 
impurities of AG. Instead, it might be scattering on ex- 
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citations of the Kondo system of interacting local mo- 
ments |8] . The situation is far from being transparent as 
is shown by transport measurements with magnetic and 
non-magnetic substitutions Alternatively, the gap- 
lessness might be caused by only a part of the Fermi sur- 
face being fully gapped in multiband scenarios (TT 
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Since our technique only provides information on the 
change AA(T), we compare our data with those of Ref. @ 
where A(0) was estimated from the surface impedance 
data. By shifting our AA(T) (circles in Fig. [I]) we obtain 
one-to-one correspondence with the data of Ref. @] shown 
in Fig.[T]by a grey band (the width is artificial - otherwise 
the data just collapse on top of each other). The solid 
line is the fit to Eq. |T]) in the full temperature interval. 
The only fitting parameter. A(0) = 358 nm, is larger than 
281 nm obtained in Ref. [6| from the low T part of the 
data (shown in the inset). The fit is of a high precision 
reflected in the value of the "coefficient of determination" 
R 2 = 0.99991 (for the perfect fit R 2 = 1). 

Inset to Fig. [I] shows comparison of the data in the 
low-temperature domain with the form suggested in 
AA = A - A(0) oc T 3 / 2 , and with the AA oc T 2 , expected 
from our model (displaced vertically by 20 nm for clar- 
ity). Straight lines show clearly that the T 2 dependence 
describes the data better. 

Long experience of dealing with pair-breaking effects 
reveals that all of them are described by the AG formal- 
ism, provided the dimensionless pair-breaking parame- 
ters are properly defined for each particular case [Hj]. 
We therefore choose the AG scheme and characterize the 
scattering by two parameters 



p = h/2irT c r , p m = hjl-nTcTn 



(2) 



where T c is the critical temperature (not to confuse with 
hypothetical T c q of the material free of scattering), 1/r 
and l/r m are the transport and the pair- breaking scat- 
tering rates. As mentioned, AG find A~ 2 cx 1 — t 2 for 
sufficiently strong pair-breaking in dirty isotropic s-wave 
materials withp 3> p m . The material of interest here is a 
clean d-wave [3j; we show below that this T dependence 
holds in this situation as well. 

Within microscopic theory, the penetration of weak 
magnetic fields into superconductors is evaluated by first 
solving for the unperturbed zero-field state and then 
treating effects of small fields as perturbations. Perhaps, 
the simplest for this task is the Eilenberger quasiclassi- 
cal approach [3]. Formally, it consists of equations for 
Eilenberger functions f(r,v,uj), /+, and g which orig- 
inate from Gor'kov's Green's functions integrated over 
the energy near the Fermi surface to exclude fast spatial 



oscillations on the scale 1/kp: 

vUf = ^ - 2ojf + ±- (/> - -L (g) , (3) 
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Here, v is the Fermi velocity, II = V + 2iriA/4> Q , tp is 
the flux quantum. A is the order parameter which may 
depend on the position fcj? at the Fermi surface (or on 
v) in cases other than the isotropic s-wave; for a d-wave 
material with a simple cylindrical Fermi surface the order 
parameter can be written as 



A = *(r,T)0, ^ = ^008 2^, 



(6) 



where ip is the azimuthal angle on the Fermi cylinder and 
f2 is normalized to have (f2 2 ) = 1. 

Further, ui are Matsubara frequencies defined by hio = 
irT(2n + 1) with an integer n, (...) denote averages over 
the Fermi surface, and 
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(7) 



The system ©-© should be complemented with the 
self-consistency equation for the order parameter and 
with an expression for the current density. For the d- 
wave symmetry, both magnetic and non-magnetic scat- 
tering suppress the critical temperature [la ]. The self- 
consistency equation in the form taking this into account, 
reads: 
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where p + = p + p m . Finally, the current density expres- 
sion completes the Eilenberger system: 



j = -47r|e|A r (0)T Im^ (vg 



LU>0 



(9) 



N(Q) is the total density of states at the Fermi level per 
one spin. 

Calculation of A(T;r, r m ) for arbitrary r's is difficult 
analytically. However, for a strong T c suppression, the 
problem is manageable. We begin with uniform zero- 
field state for which p + p m is close to the critical value 
where T c — > 0; in this state / <C 1 and g w 1 — f 2 /2 in the 
whole temperature range [7[. One can look for solutions 
of Eq. ([3]) as / = /i + /2 with f 2 <C /i . In the lowest 
approximation Eq. ^ yields: 



2A/h-2tof 1 + (f 1 )/T- -h/r^ 



0. 



(10) 



Since (A) = for the d-wave materials, averaging of this 
equation over the Fermi surface gives (fx) = as well. 
Taking this into account we have 



fi = A//^ 



1/2T 1 



(11) 
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The next approximation yields: 
A / (A 2 ) 



h = 



2/i 3 w+ 3 V 2r+w+ 



- A' 



(12) 



Making use of A = m, (0 2 ) = 1, (il 4 ) =3/4 for a 
Fermi cylinder, we evaluate: 
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Substitute this in Eq. ((8]), express the sums in terms of di- 
gamma functions, and utilize the asymptotic expansion 
ip{z + 1/2) w lnz + l/24z 2 for z = p+/2t > 1: 



■lni = tp 

3fr 2 
32tt 2 T 2 

; - hit + 



£1 
2i 

t 2 -l 
6,0+ 2 



P ' , 1 \ P .iMi ( P 



2i 2 J 9t \ 2t 

^2 



Hence, we obtain: 



48tt 2 T c V 



8tt 2 (T 2 - T 2 ) . 



(14) 
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This differs from the result for isotropic s-wave supercon- 
ductors with magnetic impurities of nearly critical con- 
centration, by a 4 times larger pre- factor [2|. The ratio 



A ma:c (0)/T c = 4^ 



(16) 



is considerably larger than 1.764 of the BCS theory for 
isotropic s-wave materials and 2.14 for the clean isotropic 
d-wave case. This ratio estimated from the tunneling 
data on CeCoIn 5 is about 10 fl6( ] . 

It is worth noting that even without magnetic scat- 
terers, the superconductivity in d-wave materials be- 
comes gapless in the domain of interest here with 
p + 1. To see this, examine the density of states 
N(e) = N(0)Reg(huj -> ie) using 5 = 1- / 2 /2 = 
1- (A/^+) 2 /2: 
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(17) 



Hence, at zero energy, N(e) has a non-zero minimum (i.e., 
the superconductivity is gapless) , whereas the maximum 
of N(e) is reached at e m — H\/3/2t + . 

Weak supercurrents and fields leave the order param- 
eter modulus unchanged, but cause the condensate to 
acquire an overall phase 9{r). We therefore look for per- 
turbed solutions in the form: 



A = A e i9 , / = (./o + /i)e ie , 
f + = {.h + .ft)e~ l \ 9 = 9 +9i, 



(18) 



where the subscript 1 denotes small corrections. In the 
London limit, the only coordinate dependence is that of 
the phase 9, i.e., fi,gi can be taken as r independent. 



The Eilenberger equations ©-([H]) provide the correc- 
tions among which we need only g\\ 
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Here P = V# + 2nA/(po = 2na/(f>() with the "gauge 
invariant vector potential" a, and 



A = A + H{f) /2r~ , Cj = uj+(g)/2i 



(20) 



In the case of interest, fo sa A/huj + -C 1, and the de- 
nominator in Eq. (fT9|) is taken in the lowest order. We 
now substitute go + g\ in the current ((9]) and compare 
the result with Anji/c = — (A 2 )^ 1 afe to obtain: 
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(21) 



Using the data of Fig.Q]with estimates for N(0) and v 
taken from we obtain p + « 8. Estimates of Q yield 
p + w 5. Hence, the statement that CeCoIns is gapless is 
not only in excellent agreement with the T dependence 
of A, but the scattering parameter p + is sufficiently large 
which is needed for our model to hold. 

There are quite a few reports of H c2 (T) for CeCoIn 5 , 
see, e.g., 17]. Our data are shown in Fig. 2. For a 
strong pair-breaking model with fixed scattering param- 
eters H C 2 cx (1 — t 2 ) [t[ that is clearly different from the 
experimental behavior. In a strong paramagnet such as 
CeCoIns, the magnetic scattering rate (and p m ) may it- 
self depend on the applied field, making our model with a 
constant p m inapplicable per se along the whole H e 2(T) 
curve. However, near T c where H C 2 — > 0, we expect the 
model to hold. We, therefore, evaluate only H c2 (T) near 
T c ; in other words, we derive the Ginzburg-Landau equa- 
tion containing the coherence length £ for the gapless case 
with a strong pair-breaking. 

Near T c , we look for a solution / of Eq. Q as an ex- 
pansion in two small parameters 



huj+ 



A_ 

TV. 



St 112 , and 



t>nA 

hw+ 2 



A 



St 



where 5t = 1 — T/T c ; smallness of the second parameter 
comes from the "slow variation" requirement. We obtain 
after simple algebra: 
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hw+ 2huj+ 2 

Further, the self-consistency Eq. which now reads as 
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should be taken into account. To evaluate (f2/), we sub- 
stitute .g = l- ff + /2 in Eq. ©, multiply it by 
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and take average over the Fermi surface: 
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Writing the last two terms, we can take / in the lowest 
order; also we make use of A = ^>fl, (il 2 ) = 1, and of 
(SI 4 ) =3/4. On the left-hand side we have: 
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summation over the repeated indices is implied. 

We now sum up Eq. (|24|) over u> and take (^/) from 
the self-consistency Eq. (|23[) . All sums obtained are ex- 
pressed in terms of di-gamma functions, for which the 
large argument asymptotics can be used. We are inter- 
ested in the field along the c crystal axis, whereas the 
plane ab can be taken as isotropic. After a straightfor- 
ward algebra one obtains the GL equation in the standard 
form 



with 
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(27) 



We note that \Pg is the zero-field order parameter of 
Eq. (|15p obtained there by a different method. 

Although the scattering parameters do not enter ex- 
plicitly the coherence length, they affect £ through 
T c (p + ). Equation (|2"7} shows that the scattering, what- 
ever it is, enhances £. In other words, in the gapless d- 
wave regime, the effect of scattering upon the coherence 
length is opposite to the familiar s-wave situation where 
£ is suppressed by scattering. When the scattering ap- 
proaches the critical value for which T c —> 0, £ diverges. 
The slope of the upper critical field at T — > T c follows: 



167T</> T c fc| 
3^2 



(28) 



where the temperature is given in Kelvins. Hence, the 
slope decreases when scattering is intensified. This un- 
usual behavior is confirmed by the data on H C 2 for low 
concentration La doping of Cei-aXa^ Coins |5j. More- 
over, we note that the Fermi velocity can now be ex- 
pressed in terms of the measured slope and T c . This 
yields v = 1.1 x 10 6 cm/s, the value between v = 2 x 
10 6 cm/s obtained using the heat capacity and resistivity 
data in Ref. [a| and v = 0.77 x 10 6 cm/s as estimated in 
Ref. [jj. 

To conclude, we show that the data on the T depen- 
dence of the penetration depth in the full temperature 




FIG. 2: (Color online) The upper critical field H c2 (T) \\ c. 
Inset shows the domain near T c with an estimate of the slope. 



range and on the upper critical field near T c strongly sup- 
port the notion that the superconductivity in CeCoIns is 
gapless. 
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